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ABSTRACT: A realistic grand unified model has never been constructed in the literature 
due to three major difficulties: the seesaw mechanism without spoiling gauge coupling 
unification, the doublet-triplet splitting and the proton decay suppression. We propose a 
renormalizable supersymmetric SO(10) model with all these difficulties solved naturally by 
imposing an extra discrete symmetry. 
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1 Introduction 

The Supersymmetric (SUSY) Grand Unified Theory (GUT) is a very important direc¬ 
tion for the physics beyond the Standard Model (SM). Among all the SUSY GUT mod¬ 
els, the renormalizble SUSY SO(IO) modelsjl, 2] are very interesting since they are very 
predictive[3-7]. All the fermions of every generation are contained in a spinor representa¬ 
tion 16 (-0), which includes also a right-handed neutrino. By coupling the ijj’s to the Higgs 
superfields in 10,126,120, correct fermion masses and mixing can be generated while the 
right-handed neutrinos are given masses by the Vacuum Expectation Value (VEV) from 
the SM singlet of 126 so that low energy neutrino oscillations are explained through the 
seesaw mechanism[8-16]. 

There are three major difficulties in the SUSY SO(IO) models to be solved. First, there 
is a conflict between gauge coupling unification, which requires a one-step breaking of the 
SO(IO) group into the SM gauge group at a GUT scale A-gut ~ 2 x 10 16 GeV[5], and the 
seesaw mechanism which requires the existence of a seesaw scale A seesaw ~ 10 ^A gut- Such 
a lower intermediate scale not only breaks the unification, but also lowers the masses of 
some color triplet masses which creates more danger proton decay [5]. Second, proton decay 
is not suppressed to satisfy the present experimental limit. The dominant mechanism for 
proton decay in SUSY GUTs is through the dimension-five operators mediated by the color 
triplet higgsinos coupling to the matter superfields. These triplets have masses at A gut 
which are not large enough. Third, there are a pair of weak doublets which are responsible 
to the electro-weak symmetry breaking in the Minimal SUSY SM (MSSM) and must have 
weak scale masses ~ 10 2 ~ 3 GeV. The large doublet-triplet splitting in masses, which is 
more serious than the splitting in the masses between neutrinos and charged leptons of 
order 10~ 6 ~ 10 -9 to be explained through the seesaw mechanism, needs to be naturally 
realized, although the non-renormalization theorem of SUSY models can be used to stabilize 
this splitting when it is generated. In the literature, the Dimopoulos-Wilczek (DW)[17- 
29] mechanism of missing VEV has been widely used in non-renormalizble models. The 
doublet-triplet splitting in 10 is realized by coupling two different 10 to 45 whose VEV is 
only in the (15,1,1) but not in the (1,1,3) direction, where the numbers in the brackets are 
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representations under the SU{4)c x SU(2) R x SU(2) r subgroup of SO(IO). Since in 10 the 
doublets are in ( 1 , 2 , 2 ) and the triplets are in ( 6 , 1 , 1 ), all the doublets are massless while their 
triplet partners are massive if the DW mechanism is applied. There is also the Compliment 
to the DW (CDW) mechanism with another 45 whose VEV is in the (1,1,3) direction only, 
coupling it with the second and a third massive 10 give masses to the second pair of doublets 
while forbidding proton decay mediated by the first 10. However, neither the DW nor the 
CDW mechanism can be applied directly in the renormalizable models in the presence of 
120 which contains also a pair of doublets in (15,2,2) whose existence invalidates the DW 
mechanism. Note that without 120 the supersymmetric SO(IO) scenario is not consistent 
with data [30]. 

In the present work we are aiming at building a fully realistic model of renormalizable 
SUSY SO(IO) solving all the above difficulties, following the very recent progresses on 
solving these difficulties by extending the Higgs sector. To naturally realize the model, an 
extra symmetry is enforced which is sufficient in constructing the required superpotential. 


2 Review on the very recent progresses 

It has been realized in [31] that to generate masses for the right-handed neutrinos, it is a 
VEV of order A seesaw which is needed, instead of a symmetry breaking scale above which 
new particles emerge breaking gauge coupling unification. A seesaw can be easily generated 
by introducing two 126s (Ai 2 ) and their conjugates 126s (A 12 ), with only Ai couples 
with the matter superfields. Introducing the superpotential 

( m A 12 + k®)AiA 2 + (mA 21 + $)A 2 Ai + QA 2 A 2 , (2.1) 

where 4> is a 210 and Q is a singlet, and we have suppressed dimensionless couplings except 
k which is the ratio of the first two trilinear couplings. Maintaining SUSY requires both 
the D- and F-flatness conditions which are 


\vir \ 2 + \v 2R \ 2 = \v 1R \ 2 + \v 2R f 


( 2 . 2 ) 


and 


0 = 


0 = 


VlR V 2R 


0 

mA 21 + $c 


mAi 2 + k<S>c 

Q 


0 m Al2 + k $0 \ / vir 
mA 21 + $0 Q ) \v 2 r 


(2.3) 


respectively, where 4 > 0 = $i(l, 1 , 1)7575 + $2(15,1, 1 ) 7^75 + $3(15,1, 3)^ is a combina¬ 
tion of the three VEVs of 4> responsible for SO(IO) breaking, and V(i i 2 ).RjII(i, 2 ) are VEVs 
of A( 1j2 ), A( 1j2 ) responsible for U(1)b~l breaking. For (2.3) to have nontrivial solutions, 
the determinant of the 2x2 matrix needs be zero. Physically, since this matrix is also the 
mass matrix of the SM singlets whose VEVs (v lR and ihr) break U{1)b-Li it must have 
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one unique pair of Goldstone modes. Then, choosing the (2,1) element of the 2x2 matrix 
to be zero, the solutions to (2.3) are 


vir = v 2 R 


Q 


km a 21 — mA, 


(2.4) 


and 


If we take 


v 2 r = 0. 


Q ~ 10 2 A gut, 


(2.5) 


( 2 . 6 ) 


the seesaw scale VEV of the same order is generated for v\ r from (2.4), and v 2 r ~ v\r 

follows (2.2). Although this VEV is put in by hand, it may be linked with the (reduced) 
A 2 

Planck scale as Q ~ Mp l VT k [32] by the Green-Schwarz mechanism[33-36], if we introduce 
a third pair of 126 — 126 which couple with a SO(IO) singlet which has an anomalous U(l) 
charge and gets a Planck scale VEV. 

This small value Q is also important in suppressing proton decay[32, 37]. We also take 
(2.1) as an example. There are three pairs of color triplets in each of Aj + Aj, and one 
more pair from <h. The triplet mass matrix is 


0 

00 

X 

00 

(AgC/t)3x4 

(Ag(7t)4x3 

C*4x4 


(2.7) 


where the triplets from <J> are put in the 4th column and the 4th row, and 


C*4x4 — 


/ M^t Vir Vir 0 \ 

vTr 0 0 Q 

vTr Q 0 0 
\ 0 0 Q 0 ) 


( 2 . 8 ) 


has three small eigenvalues of order v\r ~ Q following (2.4). Then, after integrating out the 

triplets in <h, A 2 , A 2 , we can get the effective triplet mass matrix whose three eigenvalues 

A 2 

are now all of the order W T , which are by a factor 0(100) larger than A gut so that proton 
decay is suppressed. This is an inverse analogue to the seesaw mechanism, and it proposes 
a relation between the seesaw scale mass with the suppression of proton decay. A practical 
model is presented for the Higgs sector with 10,126,126,210 in [37], and including 120 
has been done in [32]. However, in [32, 37] the MSSM doublets are given by fine-tuning the 
doublet mass matrix. 

The successful doublet-triplet splitting through the DW mechanism in the renormaliz- 
able models has been realized in [38], where the first 10 which couples to matter fields does 
not couple directly to a second 10 through 45. Instead, a filter sector is proposed since any 
singlet does not couple 10 to 120. The relevant superpotential is 

PH\h + m^hh + AhH 2 + — Mr 2 H 2 , (2-9) 
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where P is a singlet, h and h are all 10s and A is a 45 whose VEV is in the (15,1,1) direction 
only. Following a previous observational] that in a model with 10,126,126,210, if one 
takes the superpotential as 

4> (Hx + A) A + M a A A (2.10) 

plus terms containing only d>, the mass matrix for the doublets is 


0 

to 

X 

00 

^2x1 

to 

to 

X 

CO 

0 

to 

X 


( 2 . 11 ) 


for the bases as (77“, A**, A u ) in the columns and (Hf, A d , 4> d , A d ) in the rows. The 
first three columns are not independent with a combination of them gives a massless H u , 
while the first two rows make a massless Hj. The key point is the absence of H i<hA which, 
if present, gives a nonzero value to the (1,3) matrix element. Similarly there are also a 
pair of massless triplets. After applying the DW mechanism, the triplets become massive 
and the doublets remain massless. 120 can be included if a pair of them are used, and 
the filter sector (2.9) is needed[38]. In [38] a large P is needed to suppress proton decay. 
However, this makes the main components of the MSSM doublets are not from Hi, so that 
it is difficult to give the top quark a big mass. 


3 The model 

In the present work we will firstly give a pair of massless doublets and a pair of massless 
triplets analogue to the mechanism in [21, 38], then use the DW mechanism to give the 
triples masses. Since 120 is present, a filter sector is needed. In addition, we will use 
the CDW mechanism to forbid proton decay mediated by 10, so that no large VEV for 
the singlet P in the filter sector is needed. The other proton decay amplitudes mediated 
by 120,126 are suppressed by building up the triplet mass matrix analogue to (2.7) but 
extended. 

In order to suppress proton decay mediated by the 120 and 126(126), we need to 
double them as Dip and Ai 2 , Ai 2 , respectively. Besides Hi, only D\ and Ai are allowed 
to couple with the MSSM matter superfields. The superpotential for this sector is 

Wda = {k$ + mAu) AiA 2 + (<!> + mA 21 ) A 2 Ai + QA 2 A 2 

+ d>Di (A 2 + A 2 ) + <f> (Ai + Ai) D 2 (3-1) 

+ QHi (D 2 + A 2 + A 2 ) + (m£> + $) D 1 D 2 + QDl 

The D- and F-flatness conditions for v^ip)R,V(i,2) are the same as (2.2,2.3) with the seesaw 
VEV for vrr in (2.4) if Q ~ 10” 2 Agut- 

At first sight, the simultaneous existence of both <f>(ifi + Di + Ai)A 2 and $(H 1 + D 1 + 
Ai)A 2 might invalidate the observation following ( 2 . 10 ) that the absence of one of these 
couplings is the key point of generating massless doublet and triplet pairs. The subtlety is 
that V 2 R = 0 given in (2.5) as a consequence of SUSY, which eliminates the crossing entries 
between 4> and Hi + Di + Ai + Ai if they are proportional to V 2 R in the mass matrices for 
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the doublets and the triplets, while those crossing entries exist if they are proportional to 
V 2 R. The mass matrix for the doublets is 


mE a = 


06x5 

Aqx5 

^4x5 

C4x5 


(3.2) 


where the columns are (i7“, D™, D'™, A“; 4>“; A 2 , A?j, D%, D'] 1 ), while the rows are simi¬ 

lar. In (3.2), the 6th row correponds to and the first 5 entries in this row are proportional 
to V 2 R which is zero, according to (2.5). Then the massless eigenstates are 


h° u = y, _ a * xu ’ H d= E_ a x xd • ( 3 - 3 ) 

X^uDuDxAuAi X=H 1 ,D 1 ,D 1 ,AiAiA 


Note that the absence of B — L violating component <f>“ in H® suggests that there is no 
type-II seesaw contribution to the neutrino masses, as was discussed in [39]. The relation of 
large atmospheric mixing and the small quark 2-3 mixing based on the type-II seesaw[40, 41] 
is absence in the present model. 

Comparing to the doublets, there are two more pairs of triplets from Ai 2 + Ai 2 - The 
mass matrix for the triplets is 


mE a = 


07x6 

^7x6 

^5x6 

^5x6 


(3.4) 


where the columns are 

(Hi , Dj, D'E, A^, A'E, A^; A^, A^, A^, Dj, D'E), while the rows are similar. Again, 

there is a pair of massless triplets. We can re-write the mass matrix in (3.4) as 


mE a = 


06x6 

^6x6 

P<3X6 

CU 


(3.5) 


Note that in the lower right sub-matrix C ' ex6 , the (1,1) entry is the mass of the triplet 
from 4> which is ~ A gut, while the other entries in the first row are all proportional to 
v\r ~ A gut except one zero. All the other entries in the lower 5 rows of Cg x 6 are, besides 
the zeros, either proportional to v\r or Q which are of 10 _2 A gut- Consequently, there 
are 5 small eigenvalues in C' 6x6 which are not enough to generate 6 large effective triplet 
masses. The 6th large effective triplet mass will be generated after the application of the 
CDW mechanism. 

To make the massless triplet pair massive while keeping the massless doublets, we need 
to apply the DW mechanism. In realizing both the DW and the CDW mechanisms, we find 
three 45s (A, A', A") and one 54 ( E) are needed with the superpotential 


W DW = PAA' + (M A , A „ + E)A'A" 


(3.6) 


containing all possible interactions of A' and A" with fields which may have large VEVs. 
As will be seen later, the singlet P can be chosen the same as that in the filter sector. 
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Labelling A\ and A2 etc. as the VEVs in the (1,1,3) and (15,1,1) directions, respectively, 
the F-flatness conditions for the VEVs of A" are 


0 = (Ma'A" + 



0 = ( Ma’A» 



(3.7) 

(3.8) 


then either A! x or A' 2 is zero. When we take the CDW solution A 2 

li 


TjfiE = 0 so that the F-flatness condition for A\ is 


0, (3.7) gives Ma'A” + 


0 = PAi, (3.9) 

(3.9) is followed by either P or A\ is zero, and we choose the later which is the DW solution. 

Accordingly, we will introduce a filter sector which include the singlet P and a pair of 
10 s (h, h ) and use the following superpotential 


Wfilter — PH\h + mhhh, + AhH2 + A H2H3 + —M}j 3 H^. 

Then the relevant mass matrix for the doublets is 


/ 

0 

a d H P 

0 

0 

0 

\ 


a u Hi P 

0 

m h 

0 

0 



0 

m h 

0 

Ai = 

0 0 



0 

0 

Ai = 

0 0 



V 

0 

0 

0 

A 

Mh 3 , 

/ 


(3.10) 


(3.11) 


where the bases are (H®^, h, h, H2, H$) with H® d given in (3.3). Then we have a pair of 
massless eigenstates 


m h H* - (a u Hl P)h u 
^\ a u Hi P \ 2 + \ m h \ 2 
m h H° d - ( a d Hl P)h d 

^J\a d Hi P\ 2 + \m h \ 2 


(3.12) 


which are the weak doublets in the MSSM. For P having a VEV of order A gut, the 
components of H^ ,d in the MSSM doublet H u d are not small, so that there is no difficulty 
in giving the top quark a large mass, and a large tan f3 is needed for the small bottom 
and r masses. Different from the prediction of tan j3 ~ 1 got by fine-tuning the doublet 
mass matrix to have light eigenstates[37], generating these eigenstates through the DW 
mechanism exhibits the non-trivial aspect of the present work. 

Different from the treatment on the weak doublets aiming at giving predictions on 
the MSSM doublets, in the color triplet sector the effective triplet masses which determine 
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proton lifetime are our main concern. The triplet mass matrix following Wfut er is 


Mp lter = 


/ 0 P 0 0 0 \ 

P 0 rrih 0 0 

0 rrih 0 A 2 0 

0 0 A 2 0 A' 2 = 0 

V 0 0 0 A ' 2 = 0 Mh 3 


(3.13) 


where the bases are similar to the doublets. The effective triplet masses for H\ are generated 
which are infinity so that proton decay mediated by H\ is strictly forbidden. Together with 
those 5 large effective triplet masses following (3.5), all the proton decay amplitudes are 
suppressed. 

Alternatively, integrating out only the triplets in h, H 2 , H 3 while keeping those in h, the 
triplet mass matrix in (3.5) is enlarged by a 13th column which corresponds to the triplet 
in h and has only one nonzero entry P on the top, and by a 13th row which corresponds 
to the anti-triplet in h and has only one nonzero entry P in the left. Consequently, the 
lower-right sub-matrix C ' 6x6 in (3.5) is replaced by a new sub-matrix of 7 x 7 which has 
6 small eigenvalues including a zero. Then integrating out this sub-matrix induces 6 large 
eigenvalues including an infinity which are sufficient to suppress all proton decay amplitudes. 

Following (3.11,3.13) the determinant of Mp lUr is Aq UTi which is of the same order 
as 

A GUT Bet'(M f D ilter ), (3.14) 

where Det '(M^ lter ) = lim^o ^-Det (M^ lter + ef 4 X 4 ) is the product of all masses of the 
doublets except the MSSM doublets[5, 6 ], so there is no large threshold effect in the doublet- 
triplet sector relevant for gauge coupling unification. Possible threshold effects from other 
states can be included following [42-45], so that gauge coupling unification can be fully 
realized by adjusting the parameters of the model. 

The above results need to be protected by an extra symmetry. After very difficult 
efforts, we find that the Z 24 x Z 4 symmetry can be used. In addition, to maintaining SUSY 
at A cut, we need to introduce one singlet ( R ) replacing M A 'A" in (3-6) and a second 54 
(E r ). Also, a new 45 (A'") is introduced to kill redundant massless states whose existences 
break unification. Under the Z 24 x Z\ symmetry, the transformation properties of all the 
particles are listed in Table 1. 

The complete superpotential is 


W 


Full 


_ tt /-Full 1 ti/ Full 1 ti /-Full 1 tt /■ Full 

W DW + W SB + W DA + yy filter > 


(3.15) 


where 

W£$ = PAA' + (E + R)A'A" + A'A" A"', 

W^ u = + 4> 3 + + 4>A 2 + i rriE'E ' 2 + i m E E 2 

+ E 2 E' + E ' 3 + ^m A A 2 + E'A 2 + i m R R 2 + REE' 
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A'" 

A" 

E 

R 

A' 

P 

A 

E' 

4- 

Q 


Z 2 A 

12 

2 

12 

12 

10 

2 

12 

0 

0 

4 

-1 

z 4 

0 

-1 

0 

0 

1 

1 

2 

0 

0 

0 

0 


H x 

h 

h 

h 2 

h 3 

D\ 

Ai 

Ai 

D 2 

a 2 

a 2 

■^24 

2 

-4 

4 

8 

6 

2 

2 

2 

-2 

-2 

-2 

^4 

0 

-1 

1 

1 

2 

0 

0 

0 

0 

0 

0 


Table 1. Z 2 a x Z 4 properties of all superfields. Here ipi(i = 1,2,3) are matter superfields. 


+ + E'A"' 2 + <1>A" /2 , 

WEf = + m Al2 ) A 7 A 2 + ($ + m A21 ) A^Ai + QA 2 A 2 

+ TTh (A 2 + A 2 ) + $ (Ai + Ai) D2 + QD 2 
+ i? , (AiA 2 + AiA 2 ) 

+ TF/i (Z ) 2 + A 2 + A 2 ) + (rriD + F/ + <f>) D1D2 , 

fitter = PHih + (£' + m h ) hh + AhH 2 + A'tf 2 tf 3 

+ ^+ £)#!• 

Obviously, W^w gives the DW and CDW solutions for A and A' , respectively. Its last 
term does not contribute to any of the F-flatness conditions. Wgg 11 is the main sector 
breaking SO(IO). , whose first terms break U(l)j 3R x U(1)b-l i n t° U(l)y, generates 

a pair of massless doublets and a pair of massless triplets. The first three terms in Wj^ r 
generate a coupling H\{P)AH 2 at the same time forbidding the coupling between H 2 and 
D 1 through A, so that in both H i j2 the doublets are massless and the triplets are massive. 
The last two terms in Wj^ r gives masses to the doublets in H 2 - Accordingly there are 
some modifications without any important changes in the results. 

We have examined all the F-flatness conditions for the SM singlets without finding any 
conflict. There are two subtleties in these conditions. For the singlet VEV P, the condition 
is 

0 = A\A 1 + A 2 j4/ 2 , 

which is automatic following A\ = 0 and A 2 = 0, the DW and the CDW solutions, respec¬ 
tively. This is due to the existence of an accident Peccei-Quinn like U(l) symmetry [46]. 
Consequently, in other flatness-conditions, the VEV of P only appears in products with 
other VEVs of the GUT breaking fields, so that it is natural to take the VEV of P the same 
order as Agut , which makes P a harmless axion[47]. 

The second subtlety is that for the singlet Q, the F-flatness condition 

0 = V2RV2R, 

which is also automatic since V 2 R = 0. Q does not enter any other condition for keeping 
SUSY except (2.3), so that it can be given a value of the seesaw scale VEV A. seesaw 


























10 2 A gut, or be generated to be Mp^ T k °f the same order[32] through the Green-Schwarz 
mechanism [33-36] . 

4 Summary 

We have proposed in the present work a renormalizable SUSY SO(IO) with the following 
results. First, we use a seesaw VEV instead a seesaw scale so that SO(IO) breaks directly 
into the SM gauge group without spoiling gauge coupling unification. We find an important 
point ( V 2 R = 0) in generating massless states before apply the DW mechanism. Second, 
naturally doublet-triplet splitting is realized through the DW mechanism using a filter 
sector, and the DW and CDW mechanisms are very simply realized. Third, proton decay is 
suppressed successfully through the realization of CDW mechanism and through the special 
structure of the color-triplet mass matrix. Especially, the proton decay amplitude mediated 
by 10, which couples with the MSSM matter superfields with the largest coupling, is strictly 
forbidden. 

Although the present model is complicated, it has solved the main difficulties and 
is thus the first realistic SUSY GUT model. Large representations used in in the present 
model, as well as in all other renormalizable SUSY GUT models, bring in the result that the 
GUT gauge coupling blows up quickly above the GUT scale. This result can be explained 
if we take the picture that in the very early universe there was a phase transition of the 
GUT symmetry breaking. Consequently, without a very clear understanding on the details 
during this phase transition, the non-perturbative behavior of the GUT gauge coupling 
above the GUT scale may not be a real problem. 

We thank Xiaojia Li for many discussions. 
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